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Abstract. Various PDE models have been suggested in order to explain and predict 
the dynamics of spiral waves in excitable media. In two landmark papers, Barkley 
noticed that some of the behaviour could be explained by the inherent Euclidean 
symmetry of these models. LeBlanc and Wulff then introduced forced Euclidean 
symmetry-breaking (FESB) to the analysis, in the form of individual translational 
symmetry-breaking (TSB) perturbations and rotational symmetry-breaking (RSB) 
perturbations; in either case, it is shown that spiral anchoring is a direct consequence 
of the FESB. 

In this article, we provide a characterization of spiral anchoring when two 
perturbations, a TSB term and a RSB term, are combined, where the TSB term 
is centered at the origin and the RSB term preserves rotations by multiples of 221^ 
where j* > 1 is an integer. When j* > 1 (such as in a modified bidomain model), 
it is shown that spirals anchor at the origin, but when j* = 1 (such as in a planar 
reaction-diffusion-advection system), spirals generically anchor away from the origin. 



AMS classification scheme numbers: 34C20, 37G40, 37L10, 37N25, 92E20 
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Scientific literature has recorded the occurrences of spiral waves in excitable media, from 
the Belousov-Zhabotinsky chemical reaction to the electrical potential in cardiac tissue 
[1,2,4,11,16,17,22-24,26,32,34,36,37]. This last instance provides an inspiration for 
the study of spiral wave, as they are believed to be a precursor to several fatal types of 
ventricular tachycardia and/or fibrillation. [18,33,34]. 

After Barkley realized that the resonant growth transition from rigid rotations 
to quasi-periodic meandering and drifting could be explained solely through the 
inherent Euclidean symmetries of excitable media (which were incorporated in the 
various models) [1-3], the theory of equivariant dynamical systems (together with the 
general center manifold reduction theorem of Sandstede, Scheel and Wulff [13,27-30]) 
has been used on numerous occasions to explain and predict spiral wave dynamics 
[5-7,10,14,15,19,20,35]. 

One of the advantages of this approach is that it provides universal, model- 
independent explanations of observed dynamics and bifurcations of spiral waves 
that occur in a priori different situations; consider for instance the aforementioned 
Hopf bifurcation from rigid rotation to quasi-periodic meandering and drifting 
(observed in numerical simulations [4] and in actual chemical reactions [22]) or the 
anchoring/repelling of spiral waves on/from a site of inhomogeneity (seen in numerical 
integrations of an Oregonator system [24], in photo-sensitive chemical reactions [37] and 
in cardiac tissue [11]). Using a model-independent approach based on forced symmetry- 
breaking, it has been shown that anchoring/repelling of rotating waves generically 
takes place in systems in which the translation symmetry of §E(2) (the group of 
translation and rotation symmetries) is broken by a small perturbation [20], although 
the center of anchoring/repelling may not be an inhomogeneity when more than one 
such perturbations are present [7]. Similarly, some dynamics of spiral waves observed in 



Spiral anchoring in anisotropic media with multiple inhomogeneities 3 
anisotropic media {e.g. phase-locking and/or linear drifting of meandering spiral waves) 
have been shown to be generic consequences of rotational symmetry-breaking [19,25,26]. 

Many of the phenomena in which spiral waves are observed experimentally are 
modeled by reaction-diffusion systems 

nil 

^ = D.V'u + f{u) (1.1) 

where m is a fc-vector valued function of time and two-dimensional space, D is a matrix of 
diffusion coefficients and / : M'^ — > M'^ is a smooth reaction term. Their use as models 
is justified as Scheel has proved that systems of this form admit have time-periodic, 
rigidly rotating spiral wave solutions [31]. Such models are §E(2)— equivariant: indeed, 
(II. ip is invariant under the transformations 

u{t, x) I — > u{t, xi cos 6' — a;2 sin ^ -|- pi,Xi sin 9 + X2 cos 9 + P2), (1.2) 

where {9,pi,p2) G §^ x ~ §E(2) and x G [12,35]. But inhomogeneous media is 
not perfectly Euclidean: as a result, (11.11) does not provide an appropriate frame for the 
study of anisotropic media with inhomogeneities. 

Consider rather the following generalization of the bidomain equations describing 
the electrical properties of anisotropic cardiac tissue 

vt = (;{u + (3 - ^v), (1-3) 
V^^' + eg{a, £:)^x2a-2 = £^^(«, £)ux2X2, 

where m is a transmembrane potential, v controls the recovery of the action potential, 
is an auxiliary potential (without obvious physical interpretation), xi is the preferred 
direction in physical space in which tissue fibers align, e is a measure of that preference, 
g and h are appropriate model functions, and a, (3 and 7 are model parameters 
[5,8,21,25]. 
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If the tissue has equal anisotropy ratios (i.e. e = 0), (11.31) decouples into 

the FitzHugh-Nagumo equations for u and v, and Poisson's equation for ^ [25]. 

Ignoring the boundary, system (ll.3p is SE(2)— equivariant when e = 0, while it is only 

Z2-i-M^— equivariant otherwise (see [21, section 2.2] for details). 

Cardiac tissue is littered with inhomogeneities. In order to make the analysis more 

tractable, let us make the modeling assumption that the inhomogeneities consist of a 

finite number of independent "sources" which are localized near distinct sites (i, . . . ,(n 

in the plane (see [7] for a similar hypothesis). This situation could be modeled by the 

(slightly) perturbed bidomain equations 

n 

where /i = (/ii, . . . , //„) G M** is a small parameter and g^'^ are smooth functions, 
uniformly bounded in their variables. 

The symmetry-breaking induced by anisotropy and the presence of inhomogeneities 
at the origin then yields an equivariant structure on the semi-flow ^t,e.n of (11.41) . the 
existence and uniqueness of which is shown in [9]. In particular, $t,e,^ 

(El) is SE(2)— equivariant when (£:,/i) = 0; 

(E2) is Z2-i-M^— equivariant when e 7^ is small and /i = 0; 

(E3) preserves rotations around (but generically not translations) when £ = and 
all the /ij are zero except /Xj*, and 

(E4) is (generically) trivially equivariant when (e, /i) is small and generic. 

When £ = and n = 1, LeBlanc and Wulff have shown that rotating waves 
generically anchor at (or are repelled by) the inhomogeneity [20]. When e = and 
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n > 1, Boily, LeBlanc and Matsui have shown that anchoring remains generic in 
parameter wedges, but that the anchoring/repelhng is generically centered away from 
the inhomogeneities [7]. When fi = 0, LeBlanc has shown that rotating waves have 
Z2— spatial symmetry [19]. 

There is no reason to believe that combining anisotropy with multiple inhomo- 
geneities will change these results qualitatively: to wit, if a spiral anchors in non- 
anisotropic media, it should also do so in anisotropic media, but with Z2— spatial symme- 
try. However, the proofs of [7,19,20] cannot be pasted together to obtain this "obvious" 
result; some mathematical difficulties had to be overcome. 

The goal of this paper, then, is to provide a detailed analysis of abstract dynamical 
systems which share properties (El), (E3) and (E4), as well as some technical condi- 
tions which will be specified as we proceed, with the following slight modification: 

(E2') the semi-fiow is Zj.+M^— equivariant when e 7^ is small and /i = 0, where j* is 
some positive integer. 

As noted, for the bidomain model we have j* = 2, and for a reaction-diffusion-advection 
system, we have j* = 1 . The value of j* will play an important role in the subsequent 
analysis. The paper is organized as follows. In the second section, we derive the center 
bundle equations of the semi-fiow of an appropriate system near a hyperbolic rotating 
wave. We state and prove our main results for n = 1 in the third section: as long as 
J* > 1, spirals anchor qualitatively as they would in the absence of anisotropy. But in the 
case J* = 1, an unexpected result occurs: in a reaction-diffusion-advection system with 
one inhomogeneity, the tug of war between the advection current and the inhomogeneity 
ends in a tie, and spirals anchor at some point away from the inhomogeneity. To the 
best of our knowledge, this has yet to be observed either in simulations or in the lab. In 
the fourth section, we look at spiral anchoring in the presence of the most general form 
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of Euclidean symmetry-breaking. Finally, we perform a simple qualitative numerical 
experiment demonstrating the validity of our results. 



2. Reduction to the Center Bundle Equations 

Set 1 < J* G N. Let X be a Banach space, W C M x M" a neighborhood of the origin, 
and ^t,e,ti be a smoothly parameterized family (parameterized by (e, fi) G U) of smooth 
local semi-flows on X. 

Let §E(2) = C+SO(2) denote the group of all planar translations and rotations, 
and let 

a:§E(2) — > GL{X) (2.1) 

be a faithful and isometric representation of SE(2) in the space of bounded, invertible 
linear operators on X. For example, if X is a space of functions with planar domain, a 
typical §E(2) action (such as (11.21) in the preceding section) is given by 

(a(7)M)(x)=M(7-i(x)), 7G§E(2). 

We will parameterize SE(2) as follows: SE(2) = C x with multiplication given 

by {pi,(pi) ■ {P2,f2) = {e"^'-p2 + Pi,'fi + f2), V(pi,v3i), {p2,^2) e C X §^ For fixed 
^ G C, define the following subgroups of SE(2): 

C+Zj, = I (^s, ky^ I A; G Z, s G C I and §0(2)^ = { (e, 0) ■ (0, 9) ■ (-^, 0) | ^ G S^, 

the latter of which is isomorphic to SO(2) and represents rotations about the point ^. 
We will assume the following symmetry conditions on the family of semi-fiows. 

Hypothesis 1 There exists n distinct points ^i, . . . ,$,n in C such that if Cj denotes the 
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j^^ vector of the canonical basis in M", then Vm G X, e ^ 0, a 7^ 0, t > 0, 

*i,e,o(a(7)M) = «(7)*i,e,o(M) 7e C+Zj*, 

^tfl,ae,{a{-f)u) = a{-f)^t,o,ae,{u) 7^ §0(2)^^,, and 

$t,o,o(a(7)^) = a{^)<l>t,o,o{u), V7 G §E(2). 

Hypothesis [1] basically states that ^t,e,^l satisfies properties (El), (E2'), (E3), (E4). We 
are interested in the effects of the forced symmetry-breaking on normally hyperbolic 
rotating waves. Therefore, we will assume the following hypothesis. 

Hypothesis 2 There exists u* & X and Vt* in the Lie algebra o/SE(2) such that e*^ * 
is a rotation and ^t.o{u*) = a{e^*^)u* for all t. Moreover, the set { A G C | |A| > 1 } 
is a spectral set for the linearization a(e~^ )D^ifl{u*) with projection such that the 
generalized eigenspace range(P*) is three dimensional. 

In order to simplify the analysis, we only consider one-armed spiral waves, i.e. the 
isotropy subgroup of u* in hypothesis [2] is trivial. It should be noted that hypotheses 1 
and 2 hold for a large variety of spirals (such as decaying spirals), but not for all spirals 
(including Archimedean spirals) [31] li 

From now on, we assume ^t.£,^l is a semi-flow that satisfy both hypotheses, as well 
as all other hypotheses required in order for the center manifold theorems of [27-30] to 
hold. Then, for (5, /i) near the origin in R x M", the essential dynamics of the semi-flow 
^t,e,^l near the rotating wave reduces to the following ordinary differential equations on 
the bundle C x (see [5] for more details): 

p = e'^ [1/ + JP(p,p,(/^,£,/i)] 

(2.2) 

if = uj,,, + .F{p,p,ip,e,ij), 



X Even when hypothesis 1 fails, finite-dimensional center-bundle equations which share the symmetries 
of the underlying abstract dynamical systems have a definite predictive value [1.2. 19,20,22]. 
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where z/ is a complex constant, 7^ c<j,„t is a real constant, JP(p, p, 0, 0) = and 
J'^(p, p, (/p, 0, 0) = 0. Furthermore, the functions and J'^ are smooth and uniformly 
bounded in p. If {e,fi) is near the origin, we can re-scale time along orbits of fl2.2p . 
perform a simple computation and apply Taylor's theorem to get the following. 

Proposition 2.1 The symmetry conditions in hypothesis U\ imply that the equations 
l\2.^ have the general form 



eG{^it),e) + J2 1'jH.iip - Qe-'^^'\ {p - 



(2.3) 



where, without loss of generality, ipit) = t, v & C, fi = (/ii, . . . and the functions 

G,Hj are smooth, periodic in (p and uniformly bounded in p, and G is 2tx / f— periodic 
in if for some positive integer j*]^ 

As in [7], a 27r— periodic solution p^.^ of (12.31) is called a perturbed rotating wave of (12. 3p . 
Define the average value 

\PejA = ^ jjeAt)dt. (2.5) 

If the Floquet multipliers of p^^fj, all lie within (resp. outside) the unit circle, we shall 
say that [Pe,^]A is the anchoring (resp. repelling, or unstable anchoring) center of p^^^. 



§ Strictly speaking, the center bundle equations take the form 

p = J^ + e'*Y^ >^i^kn,AP.P, ij,ik,t, A), (2.4) 

where J- is the right hand side of (|2.3|) . A = (Ao, Ai, . . . , A„) = (e, /ii, . . . , /i„) £ ]R"+^, and the functions 
Ti-j^k are suitably smooth and 27r— periodic in t. However, the analysis in what follows is independent 
of the "mixed" perturbation terms Ti.j^k, when ||A|| is small enough: in particular, it is then generically 
the case that if (|2.4p has a 27r— periodic solution, so does (|2.3p . and vice- versa, and these solutions share 
the same stability. The argument depends on the Taylor expansion of an the time— 27r map (see [7] for 
an example of a similar analysis that includes the mixed perturbation terms). 
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In the following section, we perform an analysis of the anchoring of perturbed 
rotating waves of (12. 3p for n = 1 and .^i = and parameter values {e, /x) near the origin 
in M^. As noted, we will tackle the case n > 1 in section HI 

3. Analysis of the Center Bundle Equations {n = = 0) 

Equations (12.31) represent the dynamics near a normally hyperbolic rotating wave for 
a parameterized family of semi-flows satisfying the forced-symmetry breaking 

conditions in hypothesis [H We start with a brief review of spiral anchoring in the 
two cases e = and fii = 0, which were studied in detail in [7, 19,20], and then present 
the result of combining both types of perturbations. 

3.1. The Case e = 

By writing w = pe~**+ iv, this system becomes 

w = —iw + fiiH{w,w, fii) (3.1) 
where H{w, w, /ii) = Hi{w — iv, w + iv, 0, /xi). The following theorem is proved in [20]. 

Theorem 3.1 Let a = Re(Di#(0, 0, 0)), where H is as m fTOj) . Ifay^O, then for all 
/ii 7^ small enough, i\2.3^ has a hyperbolic rotating wave 

pit) = {-IV + 0(/ii)) e'\ ipit) = t. (3.2) 

The origin [p]a = Q is an anchoring center if afii < 0; it is a repelling center if afii > 0. 

In the case where the semi-flow $t,o,/ii is generated by a system of planar reaction- 
diffusion partial differential equations, the solution (13.21) represents a wave which is 
rigidly and uniformly rotating around the origin in the plane. In the case where a/ii < 0, 
the rotating wave is locally asymptotically stable. When afii > 0, the rotating wave 
is unstable (see [24] for an experimental characterization of this phenomenon in an 
Oregonator model). 
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3.2. The Case /ii = 

Since G is 27r/j*— periodic in ip (or in t, equivalently), it can be written as the uniformly 
convergent Fourier series 



The following theorem is proved in [19]. 

Theorem 3.2 Let j* > 1, or j* = 1 and g-i{e) = 0. Then, for all e ^ sufficiently 
small, the solutions of ( 1^.31) are 211 / f— periodic in time with discrete Zj* — symmetry . 

In the case where the semi-flow ^t,e.o is generated by a system of planar reaction-diffusion 
partial differential equations, the solutions of (12.31) represents discrete rotating waves in 
the physical space. 



How do these two cases interact when they are combined? We shall see that the answer 
depends greatly on the nature of j*. 

3.3. The Case f = 1 

Let Fg : M X M2 ^ C be deflned by 



where the gm are the Fourier coefficients of G found in (13. 3p . Differentiating (13.41) yields 



Generically, on the center bundle C x the path {FG(t,e),t) represents a discrete 
rotating wave around the origin with trivial spatio-temporal symmetry. Set z = p — Fq. 
Then, via (13.51) . (12.31) is equivalent to 




(3.3) 




(3.4) 



Fcit, e) = [v + eG{t, e) - sg^^{e)\ . 



(3.5) 



z = eg^i{e) + fiie'^Hi{{z + FG{t,e))e c.c, /xi). 



(3.6) 
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Generically, g-i{e) 7^ in some neighbourhood of the origin. Thus, fl3.6l) does not 
generically have a periodic solution when /ii = and e 0. Furthermore, the analysis 
of (13. 6p reduces to that of (13.11) when e = (see section \3A] and [20]). Set e = ifii. 
Then (13. 6p transforms to 

z = fiiH^{ze~'\ze'\t,i,fii), (3.7) 

where H^{w,w,t,i, fii) = ig-i{efii) + e**ifi(w + ^^(t, e/Ui)e~**, c. c, /xi) is smooth, 
27r— periodic in t and uniformly bounded in w. 

Set ai = DiH{—iv, iv, 0, 0). Near {z, i, /xi) = (0, 0, 0), the time— 27r map P of (13. 7p 
is given by 

P{z,i, Hi) = z + 2Trfii e(7„i(0) + ai^; + 0(|zp) + O(yUi) + higher order terms . (3.8) 

The fixed points of (13. 8p correspond to 27r— periodic solutions of (13. 7p . and so to discrete 
rotating waves with trivial spatio-temporal symmetry in (12.30 . Let 

B^{z,z, i, /ii) = ig-i{0) + oi\z + 0(|2;|^) + 0(/ii) + higher order terms 

be the function inside the square brackets in (13.81) . Note that -B*(0, 0, 0, 0) = and that 
Z^i-B=i,(0, 0, 0, 0) = «!, which is generically not 0. By the implicit function theorem, there 
is a unique smooth function z{e,\i\) defined near (e, /ii) = (0,0) with z(0,0) = and 

5^(z(e,/ii),z(£,/ii),£,/ii) =0 

near 2 = 0. However 2; = is not in general a solution of = for small (e, /ii), unless 
5'_i(0) = 0. Thus the generic behaviour of discrete rotating waves in (12. 3 p is to locally 
drift away from the origin. This leads to the following theorem. 

Theorem 3.3 Let a\ and g^i{0) ^ be as in the preceding discussion, with Re(a;i) 7^ 0. 
Then there exists a wedge-shaped region near (e, /xi) = (0, 0) of the form 

W = {(£:,/ii) e : \e\ < K\ni\, K > 0, fii nearO} 
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such that for all (e, /ii) G W with e ^ 0, {\2.!^ has a unique hyperbolic discrete 
rotating wave T>1 with trivial spatio-temporal symmetry and center [^^^^Ja near, but 
generically not at, the origin. Furthermore, [I^^^Ja is a center of anchoring when 
/ii Re(ai) < 0. 

Proof: Let z{e,fii) be the unique continuous function solving the equation = 
for small parameter vectors {e,fii), as asserted above. When fii = 0, any 2; is a non- 
hyperbolic fixed point of P and so, from now on, we will assume that fii ^ 0. If that is the 
case, and if e and /ii are small enough, the eigenvalues ix'i,2(^, A^i) of DP{z{e, fii) 
satisfy 

||^i,2(^, fJ'i)f = 1 + 47r/ii Re(ai) + fiiO{e, fii) ^ 1, 

since Re{ai) 7^ 0. When /ii Re(ai) < 0, the eigenvalues lie inside the unit circle and the 
fixed point is asymptotically stable; otherwise, it is unstable. 

Let K > he such that the preceding discussion holds whenever e G {—K, K) and 
let W be as stated in the hypothesis. If (e, /ii) is such that the time— 27r map (13. 8p has 
a hyperbolic fixed point 2(5, /ii) near 0, then (12.31) has a hyperbolic 27r— periodic orbit 

(t) centered near the origin. 

Let 7^ £ G {—K,K) and set e = ifii. Then {e,fii) G W, as \e\ = < K\fii\, 

and Zi^^^{t) is a 27r— periodic orbit for the pair (e, /ii), which we denote by ^^^^^(t). Since 
p = z — Fg, (12. 3p has a unique perturbed rotating wave I^^^^^, with 

K,,]a = l^e,,M - FGit,e)) dt = [2,,^Ja. 

But [2£,/xi]a = 0{e,fii) as (e, /ii) and so 'Dl^^_^ is a discrete rotating wave with 
trivial spatio-temporal symmetry and center near (but generically not at) the origin. 
The conclusion about the stability of T>1 follows directly from the hyperbolicity of the 
eigenvalues of DP. □ 
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When the parameter values stray outside of W, all that can generically be said with 
certainty is that solutions of (12. 3p locally drift away from the origin, and so anchoring 
cannot be centered there. After drifting, the spiral may very well get anchored at some 
point far from the origin, depending on the global nature of the function Hi in (12.31) . 

The analysis of the situation near the e— axis involves fixed points of (13.81) at oo, 
which correspond to traveling waves in (12. 3p [19,20]. Such an analysis lies outside the 
scope of this paper; the situation will be investigated at a later date. 

3.4. The Case f >1 

Let C^(C) and C]^(C) be the spaces of continuous and continuously differentiable 
functions from M to C, respectively, and ^j'^^"' be the space of 27r/j*— periodic functions 
of the variable t. 



Then Cl,. = {f:fe ^jf/^* n ^^(C)} and Cl,, = {f : f E ^f^^* H (C)} are 



Banach spaces when endowed with the respective norms 

I |m| |o = sup{|'u(a;) I : X G [0, 27r/j*]} and | |m| |i = | |m| |o + | |o, 

and the linear operator y : C*27r/j* ~^ ^^t^/j* defined by y{u) = in + u' is bounded, 
invertible and has bounded inverse. The nonlinear operator Tic '■ C^^^^:, x ^ 
given by 



where the gm are the Fourier coefficients of G found in (13. 3p . will play an important part 
in what follows. Note that Hg(0,0,0) = and Di7^g(0, 0, 0) = z ^ 0. Thus, by the 
implicit function theorem, there is a neighbourhood M of the origin in and a unique 
smooth function f/ : ^ ^2iT/f satisfying f/(0,0) = and HciU^e, fii),e, ^i) = for 
all (£,/ii) e N. 



Hg{u, e, fii) = y{u) - fiiHi lu-iv + e ^ 




(3.9) 
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Let Fg : M X ^ C be defined by 



FG(t,£,/ii) = 



(3.10) 



Then 



and 



y{U{e,^i^)){t)= ^iiHi (FG(t,£,/ii)e-'*,c.c.,/ii) 



Fait, e, /ii) = e** [i; + ^^(t, e, /xi) + 3^ (f/ (e, /ii)) 



(3.11) 



(3.12) 



On the center bundle CxGS*^, the path {Fait, e, fii),t) (generically) represents a discrete 
rotating wave around the origin with Zj. —spatio-temporal symmetry. Set z = p — Fq- 
Then, using fl3lll) and fl312D . flOD rewrites as 

i = /iie^* [i7i((2 + Fait, e, /ii))e-^*, c.c, /ii) - H^i^Fcit, e, i^i)e-'\ c.c, /ii))] , (3.13) 



which reduces to 



where 



z = fiiC H{ze , ze , t, e, /xi), 



(3.14) 



H{w,w,t,e,iJ.i) = Hi{w + FG{t,e,ni)e '*,c.c.,/ii) - iJi(FG(t, e, /ii)e **,c.c.,/ii) 

is 27r/j*— periodic in t. 

Near (2;, e, /xi) = (0, 0, 0), the time— 27r map P of (13.141) is given by 



P{z,e,fj.i) = z + 27i^i {ai + 0{e, fii))z + 0{e, ^ii)z + 0{\z\'^) 



(3.15) 



where ai is as defined in section 13.31 The fixed points of fl3.15p correspond to 
27r— periodic solutions of (13.141) . and so to discrete rotating waves with Z^*— spatio- 
temporal symmetry in (12. 3p . Let 

B*{z,z,e, ^i) = {ai + 0{e,^i))z + 0{e,^i)z + 0(|2:|^) 
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be the function inside the square brackets in (13.1 5p . Note that -B*(0, 0,0,0) = and 
that DiB*{0, 0, 0, 0) = ai, which is generically not 0. By the imphcit function theorem, 
there is a unique smooth function z{e,fii) defined near (e, yUi) = (0,0) with ^(0,0) = 
and 

E*(z(e,/ii),2(e,/ii),e,/ii) =0 

near z = 0. But z = is always a solution of B* = 0; hence z{e, fii) = for all small 
enough (e,/ii). This leads to the following theorem. 

Theorem 3.4 Let ai be as in the preceding discussion, with Re(ai) ^ 0. // {B,fii) 
is in a small deleted neighbourhood W^* of the origin, i\2.3\) has a unique hyperbolic 
discrete rotating wave T>1^_^ with Ijf — spatio-temporal symmetry and center [D^ ^J^i = 0. 
Furthermore, ^^^^^^ is anchored at the origin i/Re(ai)/xi < 0. 

4. Spiral Anchoring With General FESB 

The case n > 1 combines the results from the preceding section and from section 3.2 in 
[7] ; it provides a synopsis of spiral anchoring under the most general form of Euclidean 
symmetry-breaking. The results depend yet again on the nature of f. The proofs follow 
the general lines of previous work and are thus omitted, in order to avoid unnecessary 
repetitions. 

Theorem 4.1 Let j* = 1 and k e {l,...,n}. Set ak = DiHk{—iv,iv,0,0). Write 
e = fiQ. Generically, there is a wedge-shaped region of the form 

Wl = {(/io, /ii, • • • , /in) e M"+^ : < IVj, fcl/ifcl, Wj^k > 0, for j k and /i^ near } 

such that for all (e, /ii, . . . , /i„,) G with e ^ 0, the center bundle equations 

i \2.3^ have a unique hyperbolic discrete rotating wave 2)^'^, with trivial spatio-temporal 
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symmetry, centered away from C,k- Furthermore, [S)^'^]a is a center of anchoring when 
/XfcRe(aA:) < 0. 

Theorem 4.2 Let j* > 1 and k E {1, . . . ,n} . Set ak = DiHk{—iv,iv, 0,0). Generical- 
ly, there is a cone-like region of the form 

W{ = {{e, Hi,..., fin) ^ IK""^^ : kl < eo, l/^jl < Wj^klf^kl, Wj^k > 0, for j k and near } 

such that for all {e, ji) G W^, , {\2.^ has a unique hyperbolic discrete rotating wave T>^*'^ 
with Xf — spatio-temporal symmetry centered at ^k- Furthermore, [S^'^jA is a center of 
anchoring when /ifcRe(afc) < 0. 

5. Numerical Simulations 

In this section, we examine a system that models excitable media with anisotropy and 
a single inhomogeneity, which give rise to a semi-flow satisfying the FESB equivariance 
conditions described by (El), (E2'), (E3) and (E4) for j* = 1. 

The computations are carried out on a two-dimensional square domair[|| with 
Neumann boundary condition, using a 5-point Laplacian and the Runge-Kutta forward 
integrator of order 2. This naive numerical approach is used instead of a more robust 
algorithm mostly because the aim of this section is to illustrate the qualitative (rather 
than quantitative) power of theorems 14.11 and 14. 2[ 

Given appropriate initial conditions, this system sustains spiral waves. A planar 
reaction-diffusion-advection system (RDAS) is a reaction-diffusion system to which an 
advection term has been added: 

ut = DAu + Aiu^^+ A2U^^ + f{u), (5.1) 

where Ai G M2(M) for i = 1, 2, and the other terms are as in (11.11) . 



More precisely, on the domain [—30, 30] with 200 grid points to a side and time-step At = 0.005. 
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Advection terms are used to model a directed flow or current through the excitable 
medium under consideration; for instance, the modified RDAS 

ut = DAu + f{u) +e Ai{e,n)u^^+ A2ie,n)u^^ + ljg{u,\\x - ^\\'^,e, fi), (5.2) 

where {e,fi) G is small, Ai,A2 : M2(M) are smooth functions and g is a. 

smooth bounded function, could model membrane potentials in a piece of cardiac tissue 
subjected to a directed current, with an inhomogeneity located at ^. 

Let $t,e.^ denote the semi-flow generated by (15. 2p : it is SE(2)— equivariant when 
(e, /i) = 0; C-i-{l}— equivariant when e 7^ is small and /i = 0; SO(2)^— equivariant 
when /i 7^ is small and e = 0, and (generically) trivially equivariant when (e, /i) is 
small and generic |5] 

However, it is not clear if theorem 13.31 can be applied to (15. 2p since it is not yet 
known if the center manifold reduction theorems of [27-30] hold in general for RDAS; 
as such, the center bundle equations need not be of the form ( 12. 2p in the presence of 
advection. 

Yet, the following example shows that the conclusion of theorem 13.31 holds for the 
following RDAS, and so that, in some sense, (12.21) may capture the essential dynamics 
near a rotating wave when the semi-flow of the PDE has the general symmetry-breaking 
properties outlined above. 

Consider the RDAS 

Mi = i (n - - v) - 3v/2sin (^) (f)2 + Au + 0.002m^.^, 

(5.3) 

Vt = <;{U + P - JV + 02), 

where <^ = 0.3, j3 = 0.6, 7 = 0.5 and </> is defined via 

0(x) = 0.12/(xi + 10, X2 - 5^3), 



^ This can be seen by slightly modifying the proof of theorem 2.1.3 in [5]. 
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Figure 1. Anchoring in tlic RDAS equations (|5.3[) . The spiral tip path is plotted in 
black. The anchored perturbed rotating wave is shown as a gray closed curve; the gray 
square indicates the location of the perturbation center. 

with f{x) = exp (—0.00086 (x^ + Xg)) . A single integration (see figure 1) clearly shows 
the spiral anchoring away from the perturbation center. 
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